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ABSTRACT 

We determine the linear amplitude of mass fluctuations in the universe, erg, 
from the abundance of massive clusters at redshifts z = 0.5 — 0.8. The evolution 
of massive clusters depends exponentially on the amplitude of mass fluctuations 
and thus provides a powerful measure of this important cosmological parameter. 
The relatively high abundance of massive clusters observed at z > 0.5, and the 
relatively slow evolution of their abundance with time, suggest a high amplitude 
of mass fluctuations: a 8 = 0.9±10% for fl m = 0.4, increasing slightly to a 8 = 0.95 
for Q m = 0.25 and a 8 = 1.0 for Q m = 0.1 (flat CDM models). We use the cluster 
abundance observed at z — 0.5 — 0.8 to derive a normalization relation from 
the high-redshift clusters, which is only weakly dependent on Q m : a 8 Q^ 4 = 
0.78 ±0.08. When combined with recent constraints from the present-day cluster 
mass function, a 8 Q°J = 0.33±0.03, we find a 8 = 0.98±0.1 and Q m = 0.17±0.05. 
Low a 8 values (< 0.7) are unlikely; they produce an order of magnitude fewer 
massive clusters than observed. 

Subject headings: cosmological parameters — galaxies: clusters: general 



1. Introduction 



The amplitude of mass fluctuations is a fundamental cosmological parameter that de- 
scribes the normalization of the linear spectrum of mass fluctuations in the early universe 
- the spectrum that seeded galaxy formation. The abundance of massive clusters depends 
exponentially on this parameter (assuming Gaussian initial fluctuations), because a high am- 
plitude of mass fluctuations forms structure rapidly at early times, while a lower amplitude 
forms structure more slowly. The most massive systems (~ 1O 15 /i _1 M ), which take the 
longest time to form and grow, did not exist at early times if the initial amplitude of mass 
fluctuations is low, but rather formed only recently. 
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The amplitude parameter, denoted a 8 when referring to the rms linear density fluc- 
tuation in spheres of radius 8/i _1 Mpc at z — 0, is not easily determined since the mass 
distribution cannot be directly observed. As a result, this parameter is not yet accurately 
known. Recent observations suggest an amplitude that ranges in value from eg ~ 0.7 to a 
'high' value of a$ ~ 0.9 — 1. While the difference in the reported values is only around 50%, 
the impact on structure formation and evolution is much larger, since the latter depends 
exponentially on erf. The low amplitude values are suggested by current observations of the 
CMB spectrum of fluctuations (Netterfield et al. 2002; Sievers et al. 2002; Bond et al. 2002; 
Ruhl et al. 2002). However, this a 8 determination is degenerate with the unknown optical 
depth at reionization: if the optical depth were underestimated, then a 8 would be as well 1 . 
Recent observations of the present-day cluster abundance as well as cosmic shear lensing 
measurements have also suggested that erg ~ 0.7 (e.g. Jarvis et al. 2003; Hamana et al. 2002; 
Seljak 2001). However, these measures provide a degenerate relation between the amplitude 
as and the mass-density parameter Q m : a 8 Q^ ~ 0.33 (Ikebe et al. 2002; Bahcall et al. 
2003; Jarvis et al. 2003; Seljak 2001). The amplitude a 8 ~ 0.7 is implied only if Vt m ~ 0.3. 
If Vt m ~ 0.2, as is suggested by some observations (e.g. Carlberg, Yee & Ellingson 1997; 
Bahcall k Fan 1998; Bahcall et al. 2000; Wilson, et al. 2001; Ikebe et al. 2002; Reiprich & 
Bohringer 2002), then the amplitude is cr 8 ~ 0.9 — 1. Early results from the Sloan Digital 
Sky Survey (SDSS) cluster data (Bahcall et al. 2003) use the shape of the observed cluster 
mass function to break the degeneracy between the parameters and find a 8 = 0.9^Q2 and 
Vt m = 0.19^0 07. Similar results have recently been obtained from the temperature function 
of a large sample of X-ray clusters (Ikebe et al. 2002; Reiprich & Bohringer 2002). Most of 
the recent cluster normalization observations, as well as cosmic shear lensing measurements 
suggest cr 8 ~ 0.9 — 1 if Vl m ~ 0.2 (Jarvis et al. 2003; Hamana et al. 2002, and the references 
above). Combining current CMB measurements with the SDSS cluster mass function yields 
intermediate values of a 8 = 0.76 ± 0.09 and Vt m = 0.26l{]$ (Melchiorri et al. 2003). 

The evolution of cluster abundance with time, especially for the most massive clusters, 
breaks the degeneracy between cr 8 and Vt m (e.g. Peebles, Daly & Juszkiewicz 1989; Eke, Cole 
& Frenk 1996; Oukbir & Blanchard 1997; Bahcall, Fan, & Cen 1997; Carlberg et al. 1997; 
Bahcall & Fan 1998; Donahue & Voit 1999; Henry 2000). This evolution depends strongly 
on a 8 , and only weakly on Q m or other parameters. The expected abundance of massive 
clusters at z ~ 0.5 — 1 differs between Gaussian models with a 8 = 0.6 and a 8 = 1 by orders- 
of-magnitude, nearly independently of other parameters (Fan, Bahcall & Cen 1997, e.g.). 
Therefore, this method provides a uniquely powerful tool in estimating the amplitude a 8 . 



See note at end of paper 
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In this paper we use the abundance of the most massive clusters observed at z ~ 0.5 — 0.8 
to place a strong limit on <t 8 . 

2. The Evolution of Cluster Abundance 

A flat tt m = 0.3 cold-dark-matter universe (LCDM) with a 8 = 0.6 predicts ~ 10 3 
fewer massive clusters (with mass ~ 1O 15 /i -1 M ) at z ks 0.8 than a universe with a 8 = 1 
(e.g. Bahcall & Fan 1998; Bode et al. 2001). We use massive clusters observed at z ps 0.5 
and z ~ 0.8 (as compiled by Bahcall & Fan 1998) to set limits on a 8 - These three clusters, 
detected originally in X-rays by the EMSS survey (Henry et al. 1992; Luppino & Gioia 1995), 
have masses larger than 8 x 10 u h~ 1 M & within a radius of 1.5/i -1 comoving Mpc. The masses 
are determined from gravitational weak lensing observations (for two of the three clusters), 
as well as from the observed temperatures (T >8 Kev) and velocity dispersions (>1200 
Km/s) of the clusters. All clusters have a measured S-Z effect (Grego 2001; Carlstrom 
et al. 2001). These clusters have been conservatively selected, with mass measurements 
available from several independent methods, all yielding consistent results (see Bahcall & 
Fan (1998) for details, including the consistency of the mass determinations from different 
methods, and the relevant abundances at z ~ 0.6 and 0.8). Since only a threshold cluster 
mass is used in the analysis below (i.e., not individual cluster masses), a conservative mass 
threshold of 8 x 1O 14 /i _1 M is used; the clusters are well above this threshold. The resulting 
abundances are n d = 1.4lJ;£ x 10" 8 /i 3 Mpc~ 3 at z = 0.5 - 0.65, and 1.4+^ x lO^Mpc^ 3 
at z — 0.65 — 0.9. The error bars represent 68% confidence level assuming Poisson statistics 
and equal likelihood for each log(n c /). The volume searched and hence the number density 
depends on the assumed cosmology; these numbers are for an LCDM model with Vt m = 0.3. 

For purposes of comparison, the abundance of massive clusters at z ~ and at z ~ 0.38 
is obtained from the temperature function observed by Ikebe et al. (2002) (for z « 0) 
and Henry (2000) (for z ps 0.38). Here we convert our threshold mass of 8 x lO 14 /i _1 M 
(within 1.5/i -1 comoving Mpc) to temperature using the mean relation between cluster mass 
(observed directly from weak gravitational lensing) and temperature: M(< lh^ 1 Mpc) = 
0.95T(kev) 10 14 /i" 1 M Q (Bahcall & Sette (2002); see also Hjorth, Oukbir & van Kampen 
(1998); Carlberg, Yee & Ellingson (1997)). These observations indicate that cluster masses 
determined from lensing are consistent on average with those derived from X-ray temper- 
atures, with an rms scatter of ~20%. For the small extrapolation to radius 1.5/i _1 Mpc 
comoving, we use the observed cluster profile on these scales: M(< R) ~ R os (Carl- 
berg, Yee & Ellingson 1997; Fischer & Tyson 1997). This allows us to determine the ob- 
served cluster abundance at z ~ and z ~ 0.38 for the relevant mass clusters. We find 
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n d = l.ll^ x l(T 7 /i 3 Mpc- 3 at z = 0.05 (Ikebe et al. 2002) and n d = 1.7±\i x 10- 8 /i 3 Mp C - 3 
at z — 0.3 — 0.5 (Henry 2000); the error bars allow for statistical uncertainty as well as the 
uncertainty in the mass threshold. These numbers assume an Q m = 0.3 LCDM cosmology. 

The observed abundance of these massive clusters as a function of redshift is presented 
in Figure 1. In the next section, we compare the data with semianalytic predictions and 
constrain the allowed range of parameters. First, we use only the high redshift (z > 0.5) 
cluster abundance (where mass thresholds are determined by multiple methods, as discussed 
above). Then we take advantage of the full evolution of the abundance, from z ~ to z ~ 0.8, 
by combining the high-redshift result with the independent z ~ cluster normalization 
relation obtained from other studies. We discuss these analyses in the following section. 

3. Comparison with Predicted Densities 

For a given cosmological model, the expected cluster mass function can be predicted 
using recent improvements to the Press-Schechter formalism. We follow the procedure out- 
lined in the Appendix of Hu & Kravtsov (2002). For a given choice of parameters, the linear 
matter power spectrum is calculated with the publicly available CMBFAST code (Seljak 
& Zaldarriaga 1996). This allows a prediction of the mass function, using a fixed mean 
overdensity as the definition of mass. This mass can then be extrapolated to 1.5 h^Mpc 
comoving, assuming an NFW density profile. The analytic predictions were compared to 
the N-body results of Bode et al. (2001), and the resulting cluster abundances were found 
to agree within 20% out to redshift z — 1. 

The first result is presented in Figure 1, where we compare the observed cluster abun- 
dance as a function of redshift with that expected from LCDM models. Here we use LCDM 
with the 'concordance' value of Q m = 0.3 (Bahcall et al. 1999, with h = 0.72 and n — 1), 
but with different amplitudes a$ ranging from 0.6 to 1.2 (bottom to top curves). Com- 
parison of data and models shows that if Q m = 0.3, then cr 8 has to be relatively high, 
(j 8 = 0.9 ± 0.1, in order to produce the observed abundance of clusters at all redshifts. Low 
normalization values of as ~ 0.7 — which would be required for Q m = 0.3 by the present-day 
cluster abundance and the cosmic shear lensing measurements discussed in §2 — produce an 
order-of-magnitude too few massive clusters as compared with observations at any redshift. 
Therefore, for f2 m =0.3, cr 8 needs to be > 0.9. 

What if Q m is not 0.3? How does this affect the allowed range of <t 8 ? In Figure 2 we 
compare the data to the cluster evolution expected for a 8 = 1.1, 0.9, and 0.7 (top to bottom 
bands), where the width of each band covers all values of VL m from 0.1 to 0.4. It can be seen 



- 5 - 



that the best fit to the data is again a 8 ^ 0.9 — 1. Low values of a 8 < 0.7 appear to be 
excluded by the data, nearly independent of fl m ; they provide significantly fewer high mass 
clusters than observed, especially at z > 0.5. This result is consistent with the independent 
SDSS cluster mass-function at low redshift which yields <j$ = 0.91q2 an d Q m — 0-19loo7 
(Bahcall et al. 2003). The independent result obtained here from the high redshift clusters 
(Figures 1-2) provides important confirmation that a 8 is indeed high (~ 0.9 — 1). 

Finally, we use the observed cluster abundance at high redshifts (z > 0.5), independent 
of the lower redshift points, to determine the best cr 8 — Q m normalization relation from high 
redshift clusters. This was done by minimizing % 2 , using the error bars given in Bahcall 
& Fan (1998) plus an additional factor of 20% to account for uncertainties in the analytic 
prediction. The allowed 68% and 95% confidence limits are presented by the solid contours 
in Figure 3. As expected, the observed abundance of these massive clusters at high redshift 
depends mostly on a 8 (§2); the dependence on Q m is very weak. (The dependence on other 
parameters - Hubble constant and spectral index - is also weak, as is discussed below). This 
fact makes the high redshift cluster abundance method a powerful tool in constraining a 8 . 
The results in Figure 3 show that for any observationally acceptable range of fl m , from ~ 0.1 
to ~ 0.4, the amplitude remains in the range of erg ~ 0.9 — 1. Values of a 8 < 0.8 are unlikely; 
they have too little power to form massive systems at z > 0.5. 

We use the z = 0.5 — 0.8 clusters to determine the best-fit cluster normalization relation 
from the high-redshift objects; we find: 



(where a 8 and fl m refer, as before, to their z = value). The linear relation [la] is within 
3% of the true best-fit a 8 over the range Q m = 0.1 — 1.0; the power-law fit [lb] is superior 
to equation [la] in the range Q, m = 0.2 — 0.7. These relations assume a Hubble constant of 
h = 0.72 and spectral index n — 1, but the results are insensitive to reasonable changes in 
h and n. We find that changing h by Ah = ±0.13, i.e., from h = 0.59 to h = 0.85, changes 
<r 8 by less than ±5% (at a given Q m ); <r 8 increases slightly with h. Similarly, changing the 
spectral index n = 1 by An = ±0.2, from n = 0.8 to n — 1.2, changes a 8 by less than 
±5% (in the same direction: slightly higher a 8 for higher n). The 1- and 2-sigma contours 
of the allowed a 8 - Q m parameter region that includes both these h and n variations are 
presented by the dotted curves in Figure 3. As expected, these conservative ranges in h 
and n broaden the allowed region. The above results include the estimated uncertainty in 
the mass threshold for the clusters (Bahcall & Fan 1998). As discussed in §2, only a lower 
threshold (selected at the lower 1-sigma level of the mass estimates) is used in this analysis, 
not the individual mass of each cluster. Reducing the cluster mass threshold further by 10% 



a 8 = 1.03 - 0.3fi m ± 10% 
a 8 n° m M = 0.78 ± 0.08 



,or 



(la) 
(lb) 
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will reduce the amplitude a 8 by ~10% (Fan, Bahcall & Cen 1997; Bahcall & Fan 1998). 

Also shown in Figure 3 are the analogous confidence contours from the cluster mass- 
function normalization at z — 0.1 — 0.2 (Bahcall et al. 2003); the best-fit relation is approx- 
imated by 

a s n° m 6 = 0.33 ± 0.03 (2) 

which has a steeper Q m dependence than the high-redshift constraint. These two independent 
constraints overlap only at low Q m . Requiring that both of the constraints of equations [1] 
and [2] be simultaneously satisfied yields 

<7 8 = 0.98 ± 0.1 (3a) 
tt m = 0.17 ±0.05 (3b) 

for the allowed 1-sigma overlap region when h—0.72 and n—1. Allowing for variations in 
h and n as discussed above yields the same central values for a 8 and Q m , but broadens the 
allowed range. Even with the very broad ranges adopted for h and n, the conservative limits 
os > 0.70 and fl m < 0.36 can be set at the 95% confidence level. 



4. Discussion and Conclusions 

We use the observed abundance of high-mass clusters of galaxies at z — 0.5 — 0.8 to 
determine the linear amplitude of mass fluctuations, a 8 . The cluster abundance depends 
exponentially on this amplitude, and only weakly on other parameters; it therefore provides 
a powerful method for measuring this important parameter. We show that the relatively high 
abundance of massive clusters observed at z > 0.5, as well as their relatively slow evolution 
with time, requires a high amplitude of mass fluctuations, cr 8 ~ 0.9 — 1. This conclusion is 
nearly independent of the exact value of Q m (in the typical range of Q m ~ 0.1 — 0.4). 

We use the observed abundance at z > 0.5 to determine a normalization relation from 
high redshift clusters. The relation depends only weakly on Q m : crgfi^ 14 = 0.78 ± 0.08; 
alternatively, a linear relation of the form o~ 8 = 1.03 — 0.3Q m (±10%) provides a similarly 
good fit to the data. These fits illustrate that a 8 > 0.8 for any Q m < 0.4. For the typical 
observationally suggested value of Vl m ~ 0.2 — 0.3, the amplitude is cr 8 = 0.95 ± 0.1. We 
emphasize that this high o 8 value indicated by the cluster abundance at high redshift is 
nearly independent of the exact value of fl m . 

We combine the high redshift constraint above with the independent normalization 
relation obtained from low redshift cluster abundance — a relation that is steeper in Vt m 



-7- 



(~ fi^; 6 ; equation [2]). The combination breaks the degeneracy between the two parameters. 
We find a 8 = 0.98 ± 0.1 and Q m = 0.17 ± 0.05 (Figure 3; equation [3]). 

The high value of cr 8 required to explain the high abundance of the most massive clusters 
at z > 0.5 is consistent with the present day cluster mass function if the mass density 
parameter is low, Vl m ~ 0.2. If f2 m =0.3, then the high redshift clusters still require a high 
amplitude (erg = 0.92 ± 0.1), since this constraint is nearly independent of Q m ; the low 
redshift cluster abundance is consistent with this value at the 2-sigma level. 

These results improve upon Bahcall & Fan (1998) by using improvements over the stan- 
dard Press-Schechter formula (which does not accurately reproduce results from cosmological 
simulations at high redshift), allowing for changes in h and n, and by using a more recent 
cluster normalization relation at low redshift. The current results yield slightly lower values 
for the cosmological parameters than the previous work (the latter suggested cr 8 = 1.2 ±0.22 
and Q m = 0.2+°;^ (68%) when using the most massive clusters) but are consistent with the 
new values within the error-bars (Bahcall & Fan 1998; Fan, Bahcall & Cen 1997). 

The excess CMB fluctuations detected on small scales by the CBI (Mason et al. 2002) 
and the BIMA (Dawson et al. 2002) experiments implies (if correctly interpreted as being 
due to the S-Z effect from distant clusters) that <t 8 = 1.04 ± 0.12 (95%; Komatsu & Seljak 
2002). This is in excellent agreement with our current conclusions. We note, however, that 
this high amplitude is inconsistent with the lower value of a 8 ~ 0.7 suggested by current 
CMB data on large scales (which is degenerate with the unknown optical depth). Future 
CMB observations should clarify this current inconsistency. If massive clusters exist with 
relatively high abundance at high redshifts, as suggested by the data used here (as well as by 
deep X-ray surveys, e.g. Rosati, Borgani & Norman 2002), then these clusters should indeed 
produce the excess S-Z fluctuations observed by the CMB data. 

The relatively high abundance of massive clusters observed at z > 0.5 provides one of 
the strongest arguments for a high amplitude of mass fluctuations, cr 8 ~ 1. 

Note added March 7, 2003: The recent CMB anisotropy spectrum released by the 
WMAP team in February 2003 nicely confirms the results presented here. The constraint 
from the CMB alone is cr 8 = 0.9 ±0.1 (Spergel et al. 2003), in full agreement with the current 
high redshift cluster constraint. 

This research was supported by the National Computational Science Alliance under 
NSF Cooperative Agreement ASC97-40300, PACI Subaward 766. 
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Fig. 1. — Evolution of cluster abundance with redshift, for clusters with mass Mi.5 jCom > 
8 x lO 14 /t _1 M (within a comoving radius of 1.5/i _1 Mpc). Dots with error bars are the 
data as described in the text. The lines are the predicted number density assuming a flat 
fl m = 0.3 cosmology; each line is labeled with the as used. 
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Fig. 2. — Same as Figure 1, but for 0.1 < Vt m < 0.4. Each band is labeled with the as used 
to predict the number density, and within each band Q m is varied from 0.1 (bottom of band) 
to 0.4 (top); h = 0.72 and n = 1. 
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Fig. 3. — Confidence contours in the Q m — as plane. Solid lines are the 68% and 95% 
confidence contours found from the high-redshift cluster abundance (contours extending to 
high Q m ) and from the low-redshift SDSS HMF clusters (extending to lower <r 8 ); these assume 
h = 0.72 and n — 1 and two degrees of freedom. The shaded region satisfies both constraints 
at the 68% level. The dotted lines are the 68% and 95% limits with four degrees of freedom 
when allowing 0.59 < h < 0.85 and 0.8 < n < 1.2. 



